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Abstract 

We formulate and study the set of coupled nonlinear differential equations which 
define a series of shape invariant potentials which repeats after a cycle of p iterations. 
These cyclic shape invariant potentials enlarge the limited reservoir of known analyt¬ 
ically solvable quantum mechanical eigenvalue problems. At large values of x, cyclic 
superpotentials are found to have a linear harmonic oscillator behavior with superposed 
oscillations consisting of several systematically varying frequencies. At the origin, cyclic 
superpotentials vanish when the period p is odd, but diverge for p even. The eigenvalue 
spectrum consists of p infinite sets of equally spaced energy levels, shifted with respect 
to each other by arbitrary energies wq, wi,... ,ujp-i. As a special application, the en¬ 
ergy spacings ujk can be identified with the periodic points generated by the logistic 
map Zk+i = ’’’ZkiX ~ ^k)- Increasing the value of r and following the bifurcation route 
to chaos corresponds to studying cyclic shape invariant potentials as the period p takes 
values 1,2,4,8,... 
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The number of analytically solvable eigenvalue problems in quantum mechanics is rather 
limited. In this regard, the concept of shape invariance in supersymmetric quantum 
mechanics |Q, has proved to be very useful because it leads immediately to exactly 
solvable potentials. In this paper, we define and study cyclic shape invariant potentials. 
The supersymmetric partner Hamiltonians correspond to a series of shape invariant po¬ 
tentials which repeats after a cycle of p iterations. Such potentials have an infinite num¬ 
ber of periodically spaced eigenvalues. More precisely, the level spacings are given by 
Wo,• • •,Wo,• • • 5 Wp_i,cuo, wi,..., (see Fig. 1). The ground state is at zero en¬ 
ergy; the next (p—1) eigenvalues are Ei = J2k=o : 1 = 0,1,..., {p—2) , and all other eigen¬ 
values are obtained by adding arbitrary multiples of the quantity flp = wo -|-wi • ■ + u}p-i. 
The general formula for the excited energy levels is 

i 

nO.p + '^uJk ; {n = 0,1,2,... ,oo ; / = 0,1,... , (p - I)} . (I) 

k=0 

Effectively, the spectrum is obtained by starting with p infinite sets of energy levels with 
equal spacing Qp and shifting these sets by successive arbitrary amounts a;o,a;i,... ,ujp-i. 
The superpotentials corresponding to cyclic shape invariant potentials turn out to 

have rather interesting behavior at large x. One finds a linear dependence Q,pxj2p corre¬ 
sponding to a harmonic oscillator of angular frequency l^p/p with superposed sinusoidal 
oscillations which become more closely spaced with increasing x. We find that these oscilla¬ 
tions at large x are given by the functional form bpj sm{qpjx‘^). For any choice of the period 
p, there are several oscillation “frequencies” qpj, with j = 1, 2,..., [p]/2, where [p] is the 
largest even integer less than p. These “frequencies” are given by qpj = Qptan{Trj/p)/2p. 
The amplitudes of the oscillations bpj depend on the energy spacings wo,a;i,.... The su¬ 
perpotentials also have unusual behavior at small x. For even values of p, one finds that 
W^'P\x) diverges like 1/x, whereas for odd values of p, IF(p^( 0) = 0. Our plan is to first 
review some basic ideas of supersymmetric quantum mechanics and the concept of shape 
invariance for a change of parameters oi = /(oq). Next, we formulate the general set of 
coupled nonlinear differential equations governing cyclic shape invariant potentials of ar¬ 
bitrary period p. We describe some real transformations / which have the property that 
f^{ao) = oq, that is the parameters ak , k = 0,1,..., {p — 1) repeat with period p. We 
study and discuss the solutions of the cyclic shape invariance conditions. The special cases 
of periods p = 3,4 are treated in some detail. Many properties generalize to arbitrary values 
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of p. Numerical solutions for W^^\x) and the corresponding potentials V^\x) illustrate 
the general results. As an application, we examine cyclic shape invariant potentials corre¬ 
sponding to energy level spacings ojk which are the periodic points generated by the logistic 
map Zk+i = rzk{'\- — z^) 0]. As the parameter r increases from 1 to 3.569946, the period 
p also increases and the energy level spacings become increasingly more chaotic. This is 
reflected in the fact that the potentials V-{x) show increasingly more oscillation frequencies 
as one follows the period doubling route toward chaos. 


General Formalism for Cyclic Shape Invariant Potentials: 

Recall that the superpotential W (x, oq) generates the supersymmetric partner potentials 

V±{x,ao) = W\x,ao)±W'{x,ao) , (2) 

where oq is a set of parameters. These partner potentials are shape invariant if they both 
have the same x-dependence upto a change of parameters ai = f{ao) and an additive 
constant R{ao). The shape invariance condition is |j|] V+(x,ao) = V-{x,ai) + R{ao) , or 
equivalently 

W\x,ao) + W'{x,ao) = W\x,ai)-W\x,ai) + R{ao) . (3) 


The property of shape invariance permits an immediate analytic determination of energy 
eigenvalues |]^, eigenfunctions ||5j and scattering matrices [^. For unbroken supersymmetry, 
the eigenstates of the potential V- (x) are: 


n—1 


4 ^ = 0 , 4 ") = ^ R{ak) , 

k=0 


4^ oc e 4^ 4 \x,ao) = 


-^ + W{x,ao) 


(n = 1,2,3,...) . 


( 4 ) 


The lowest lying eigenstate is at zero energy, characteristic of unbroken supersymmetry. 

Solutions of the shape invariance condition (^) corresponding to a translational change 
of parameters ai = /(uq) = oq + /3 , = /"(oo) = + nP are well-known |p. Also, the 

scaling change of parameters ai = /(oq) = qao , = /"'(no) = Q^clq has been investigated 

and shown to produce new classes of shape invariant potentials ||^, which included self¬ 
similar potentials Q as a special case. 
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Very recently, another type of change of parameters defined by /^(oq) = ao was studied 
in detail |10|. Here, the parameters repeat after a cycle of two iterations and the spectrum 
consists of two shifted sets of equally spaced eigenvalues. In this paper, we consider a much 
broader change of parameters defined by 


/^(®o) — 0 - 0 , 


( 5 ) 


where the parameters will repeat after a cycle of p iterations. For this case, one has ag = 
) /(®o) — — ®p+i) etc. 

For a cyclic set of parameters, we now formulate the corresponding shape invariance 
conditions @ for any period p. Dehning R{ak) = for A: = 0,1,... , (p — 1), and using 
o-p = ag, one obtains a set of p coupled nonlinear differential equations 

W‘^{x,ak)+ W'{x,ak) = W‘^{x,ak+i)-W'{x,ak+i)+u;k , k = 0,1,... ,{p - 1) . ( 6 ) 

There are p unknown superpotentials W{x,ak) which are related to each other by cyclic 
permutations of the indices 0,1, 2,... , (p — 1). We are seeking superpotentials with odd 
parity, which will yield symmetric potentials V-{x). It is convenient to denote W{x,ak) = 
Wk- The set of shape invariance conditions now has a deceptively succinct form: 

Wl + W'k = Wl^^-W'k^^+u:k , k = 0,l,...,{p-l) . ( 7 ) 

There are two important problems one needs to solve: (i) to find real transformation 
functions f{y) which satisfy eq. (I) and hence produce a cyclic set of parameters ; (ii) 
to determine the superpotentials {x) by solving the cyclic shape invariance conditions 
(|^ corresponding to an arbitrary cyclic choice of parameters. These two problems will be 
treated in turn. 


Transformations Giving Cyclic Sets of Parameters: 


We now show that a simple algorithm for generating non-trivial real functions / : R 
satisfying /^(y) = y , Vy G R is the projective(Moebius) transformation 

./X ay + h 

Jyy> =-TT ’ 

cy + d 

with specihe constraints on the parameters a, b, c, d. Let M be the real 2x2 matrix 

/ 

M = 


a b 
c d 


R 

( 8 ) 

(9) 
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For f{y) given by eq. (| 8 |), observe that f^{y) is also a projective transformation whose 
coefficients are the same as the elements of the matrix M^. Therefore, solving eq. (|^ can 
be rephrased in terms of finding a real 2x2 matrix M, such that for any integer p > 1, 
the matrix is proportional to the unit matrix 1 . Without loss of generality, we choose 
the condition = 1, and our problem is to find real 2 -dimensional representations of the 
cyclic group of order p. This group being Abelian, its 2-dimensional representations are 
reducible, i.e. there exists a non-singular matrix Q for all matrices M such that 

/ 


M = Q-^MQ = 


Ai 0 
0 A 2 

Here A^ {k = 1, 2) stands for any two of the p-th roots of unity: 


( 10 ) 


V 


/27rjfc . 
= exD - 1 


Xk = exp 


V P 


jk = 0,...,p-l . 


( 11 ) 


Since we are interested in non-trivial matrices, we will impose two additional constraints: 
Ai 7 ^ A 2 and M” 7 ^ 1 for any n < p. The range of possible solutions is constrained by the 


invariance under similarity transformations (IT) of both the trace and the determinant: 


a + d = exp 


/27rji . 


V P 


i -|- exp 


27rj2 . 

- 1 

p 


; ad — be = exp 


p7r(ji +j2) , 

V P 


Also, the reality condition for the elements of M gives 


sm 


27rji \ 

p ) 


, . ( 27rj/2 A 


= 0 ; sin 


2vr(ji + j 2 ) 


= 0 . 


( 12 ) 


(13) 


Taking into account all the above considerations we are left only with the roots of unity Ai 
and A 2 satisfying ji + j 2 = p. Therefore, the matrix Ad will be a matrix of the form ( 0 ) 
with Ai = exp = A^^. Note that the choice j = 0 gives Ad = 1 , which gives a trivial 

M. Similarly, if p is a non-prime number, say p = jn, then Ad”" = 1 for n < p, which gives 
M"' = 1 , and is not a non-trivial result for M. However, for all prime periods p, the matrix 
M has the form (^) whose elements must be real and satisfy the system of equations: 


ad — 6 c = l; a-|-d = 2 cos 


(?) 


(i = 0 ,l,...,p-l) . 


(14) 


Finally, the desired transformation function f{y) is given by eq. dH) with real parame¬ 
ters a,b,c and d satisfying the constraints ([l^). For example, a projective transformation 
corresponding to the choices p = 3,j = 1 is 

ay + b 


fiv) = 


-y{a? + a + l)/b - {a + 1) 
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and it is easy to check that f^{y) = y. 


Solutions of the Cyclic Shape Invariance Conditions: 


Having established that there exist real transformations / which produce cyclic sets of 
parameters, we now turn our attention to studying the solutions of the cyclic set of equations 
(|^. A very general result comes from adding all the equations. It gives X]fc=o ~ f^p/2, 
which on integration yields 

p-i ^ 

Wk{x) = -QpX . (15) 

k=0 


This condition is valid for all x. There is no constant of integration in eq. (15) since we 
want the superpotentials Wk to be odd functions of x. In particular, at large x, constraints 
(|^ and (15) give all the linear asymptotic behaviour Q,pxj2p, corresponding to a simple 
harmonic oscillator potential of angular frequency ^pjp. Also, at small x, in order to satisfy 
eqs. (0), all superpotentials can at most have a 1 /x singularity. 

The general analytic solution of equations (^ is not easy to obtain. However, a numerical 
solution can be obtained via the Runge-Kutta method and we will also obtain several 
analytic results at small and large values of x. For ease of presentation, we systematically 
consider small values of the cycle period p before coming to general results. 


p = 1: Here oi = /(no) = oo and all energy levels are equally spaced. As expected, the 

shape invariance condition (^) gives ~ ojqxI2 which is a simple harmonic oscillator 

superpotential. 


p = 2: This is the case of two shifted sets of equally spaced eigenvalues considered in 

detail in ref. [^]. The shape invariance conditions (0) are: 

iy2 + Wq = Wl - Wl + wo , + Wi = -Wq + uji . (16) 

(2) 

The solution is obtained by straightforward manipulations. The superpotential Hq ' = 
W^‘^\x,ao) and the corresponding potential vY {x) obtained using eq. (I) are given by 
(2), . ^ 2 X (wo-wi) ( 2 ) (wo-Wi)(3wo+Wi) Wl HaV 

1^0 W-—+ , I/- T + 


Note that for period p = 2, the superpotential Hq diverges like a^/x near the origin. 


which produces a singular inverse square potential in the transition region |11, 12, 13|. A 
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detailed discussion of the potential as well as the eigenfunctions ipn ^ obtained from eq. (|^ 
is given in ref. [^]. 


p = 3: There are now three shape invariance conditions (0): 

+ Wq = wl - W[ + UJQ , + Wl = wi -W 2 + UJ 1 , + W2 = - Wo + W2 . (18) 


These equations are not easy to solve analytically. In order to obtain a numerical solution, 
it is best to hrst solve eqs. (1^) for the derivatives Wq, W{, W 2 . One gets 


Wg — wf — IV 2 T 2 ~ ’ 

W( = W 2 — Wq + ^('^1 — a;2 + ujq ) , 

W 2 = Wq — W^ + —(u^2 ~ ^0 T • 


(19) 


In this form the Runge-Kutta method is immediately applicable, and the superpotentials 
can be determined for any choice of the energy spacings u/o,uji,u’ 2 - An illustration of the 
results is shown in Fig. 2. The overall average behavior of Wq at large x is indeed 
seen to be as expected. On this linear dependence are clear oscillations which get 

more closely spaced as x increases. They correspond to a behavior bsm{qx‘^) where the 
amplitude b and the frequency q both depend on the choice of energy spacings ti;o,u;i,a; 2 . 
In particular, when one chooses ujq = toi = uj 2 , the numerical solution shows no oscillations. 
This is reasonable since for this case one expects a simple harmonic superpotential VIqx/Q. 

One can get a number of analytic results for eqs. (^) at small x as well as large x. At 
small X, if one feeds the series expansions 


Wk = ak/x + (5kX + + 5kX^ H- 


( 20 ) 


for k = 0,1,2 into eq. (19) and equates coefficients of x, one gets 


tt7-(3) ^ / \ ^ / \ ^ / 2 2 \ 

Wq = — uji -\- u: 2 ) + —ct;i(c(;o — UJ2) + ~ ~ ‘^2) + 


( 21 ) 


(3) (3) (3) 

with Wl , W 2 obtained by cyclic permutations of the indices 0,1,2 . Wq vanishes at 

X = 0. Note that J2k=olk = J2k=o^k = 0 in agreement with eq. (p^. At large x, one 

finds oscillations with one frequency qsi = This is a special case of the result for 

arbitrary p which will be derived below. 
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p = 4: It is worth discussing period 4 briefly, since several results are different from period 
3. There are four cyclic shape invariance conditions: 


+ Wq = Wl - Wl + too , Wf + ITi = Wi -W 2 +UJ 1 , 
wi -I-W 2 = wi - IT 3 + ^2 , wi + IT's = - IT^ + ^3 . 


( 22 ) 


Unlike the period 3 case, one cannot directly solve for VUg, W[, IU 2 , IU 3 , since their coefficients 
in eqs. (^) have a zero determinant. Indeed, this situation occurs for all even periods p. 
For the case of period 4, one can proceed by noting that eqs. (22) give the constraint 
Wq — VFf + IFf — IF 3 = (wg — Wl + a ;2 — W 3 )/ 2 , which on differentiation yields 

WqWq - IUiIU( + W2W2 - IU3IT3 = 0. 


This equation, along with any three of eqs. (^), permits one to solve for VFg, ITi', IT 2 , IT 3 , 
and proceed with a numerical determination of the superpotentials by the Runge-Kutta 
method. 

Analytically, at small x, substitution of expansions of the form ( 20 ) into eqs. (^) and 
collecting powers of x gives 


w(4) _ - Wl + a;2 - W3) ^ xn4(wgW3 + WgWi +W2a;3 - W2W1) ^ 

— -tt::- 1 - 77 -;- 77 -;-^-r 


(23) 


20 , 4 X ' 4 (wg + W2)(wi + W3) 

with the other superpotentials , 11 ^ 2 ^^) ^ 3 ^^ obtained by cyclic permutations of the 

indices 0,1,2,3. Note the singular 1/x dependence, similar to period 2. This is a feature 
of all even periods. At large x, period 4 superpotentials show oscillations with a single 
frequency qn = U 4 / 8 . We will derive the general result for arbitrary p below. 


Results for Arbitrary p: The experience gained from small values of the period p can 
be used to derive a number of general analytic statements valid for any arbitrary period p. 
At small X, the superpotential IFq^^(x) has the form 


ITg^^ — — (wg — Wl + • • • + Wp_i) 




3 L 


Wi (w2 — W3 + • • • + Wp_i — Wg) + W3 (w4 — W5 + 


+ • • • + Wp _2 (Wp_i — Wg + Wl — • • • — Wp-g) + 0 (x®) ; 


lui") = 


(Wg-Wi+W2 -Wp_i) + 


2UpX 


1 - 2 X 


+ Wp_i — Wg + Wl — W2) 
{p odd) 












(24) 


UJi {iU2 + <^4 + • • • + 2 ) + ^3 + ^6 + • • • + iOp— 2 ) + • • • + UJp—^U!p—2 

(wo + 0^2 + • ■ ■ + ^p-2) (<^1 + W3 + • • • + Wp-l) 

+ O(x^) ; {p even) 

The behaviour of the other superpotentials W^\x) can be seen by cyclically permuting 
the indices — 1). Note the 1/x singularity at the origin for even periods. As 

for the case of periods p = 2,4, it can be readily verified that the coefficient of the 1/x^ 
singularity in the potential lies in the transition region —1/4 to 3/4. In this region, 
both solutions of the Schrodinger equation are square integrable at the origin B Sill. 
The infinite discontinuity in W^^\x) at the origin for even p produces a 6(x) term in the 
potential vl^\x) as described in ref. (Bl)- 

At large x, one removes the asymptotic linear behaviour by defining 


Qk(x) = Wk(x) - flpx/2p , k = 0,1,... ,p - 1 . 


(25) 


Since numerical solutions suggest sinusoidal oscillations in x^, we make a change of variables 
to n = In terms of u, the equations satisfied by Qk are obtained by substituting eq. 
(^) into (0). Tetaining only the leading powers of u yields 


Q'o + Qi — i~Qo + Qi) 


(26) 


and cyclic permutations. This set of equations can be written in a more compact form by 
putting all the superpotentials into a column vector Q and making use of the matrix Cp 
whose only non-zero elements are (Cp)i 2 = (C'p )23 = ••• = (Cp}p-i^p = (Cp}pi = 1. The 
equations now read 


n 


(Cp + l)Q'-^(Cp-l)Q = 0 . 


(27) 


The most general solution is a linear combination of eigenstates of the form Q oc exp (Art 
The eigenvalues A are determined from 


det[A(Cp + 1) - (np/2p)(Cp - 1)] = 0 . (28) 

Defining A = 2pX/Q,p, and expanding the determinant, eq. ( [2^ ) can be re-written as 

(A + l)P-(-l)^’(A-l)?’ = 0 , 

which is a polynomial equation of degree p for p odd and of degree p — 1 for p even. This 
situation occurs because eq. (|B) is more general than the standard eigenvalue equation. 
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The solution is 


1 - A f2Trj \ 

—^ = exp -^) , j = 0,±1,±2,-• • ,±[p]/2 . 

\p] is the largest even integer less than p. This gives rise to distinct oscillation frequencies 

Qpj = (f^p/2p) tan(7rjyp) , j = 1,2,..., \p]/2 . (29) 

At large x, the superpotential W^\x) has the form 

^ ^ sin^qpjx"^ + (j)^) + 0{l/x) (30) 

^ i 

where the phases are 

4>P = 4>j + 2'Kk/p . (31) 

Note that the requirement of having odd superpotentials can be satisfied by inserting an 
asymptotically constant odd function such as tanh(Q;a:) or tan“^(Q;x) with the coefficient 
bpj. This does not affect the above discussion of oscillation frequencies. Also, if one uses 
eq. (^) to get the asymptotic behaviour of wave functions, the linear term in gives the 
familiar Gaussian dependence of a harmonic oscillator and the sinusoidal correction terms 
in trip) give Fresnel integrals producing oscillations on the Gaussian. 

As an example, for period 5 , taking fls = 6.5, we predict two oscillation frequencies 
(751 = 0.472 and (752 = 2.001 from eq. (|2^) , and we have checked that these values are in 
agreement with the numerical solutions shown in Fig. 3. 

Approach to Quantum Chaos Using the Logistic Map: 

We have now established the method of generating cyclic shape invariant potentials. 
This method can be applied using any transformation f{y) which generates parameters 
which repeat with period p. For the special projective transformations described previ¬ 
ously, the initial parameter oq can be chosen arbitrarily. There are also other possible 
transformations in which cyclic parameters are obtained if one starts from specific choices 
for oq. One such interesting possibility, which we discuss in detail here, is to use a set of 
cyclic parameters given by the periodic points of the logistic map 

Zk+i = rzk{l- Zk) . (32) 

This famous transformation has been extensively studied |^. For 1 < r < 3, one has a 
single non-zero fixed point aX z = {r — l)/'r; for 3 < r < 3.449499, there is a stable 2-cycle 
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with periodic points at z = [(r + 1) ± \/(r — 3)(r + l)]/2r; for 3.449499 < r < 3.544090, 
there is a stable 4-cycle, etc. As r is increased, the period keeps doubling, giving finally the 
period 2°° at 3.569946. For r > 3.569946, one obtains all odd periods and their harmonics, 
and at r > 4, there is complete chaos. 

In our application, we will identify the energy spacings with the periodic points 
of the logistic map. For example, if we choose r = 3.20, the periodic points are 0.799 
and 0.513, and we determine the potentials which have energy spacings ujq = 0.799 and 
LVi = 0.513. The corresponding superpotential is shown in Fig 4(a). Other examples with 
larger values of r and correspondingly larger periods p = 4,8 are also shown in Figs. 4(b) 
and 4(c). Clearly, as one advances along the period doubling route to chaos, more and more 
oscillation frequencies come into play - which is certainly a plausible result. Specifically, for 
period 2”, there are 2”“^ — 1 frequencies given by eq. (^). 

Although we have so far chosen values of r which yield stable cycles of period 2"', it is 
known that there are small windows of r-values where stable cycles of periods 3,5,.. .also 
exist. The above methods immediately permit us to obtain superpotentials for these cases 
also. These issues and deeper connections with chaotically spaced energy levels will be 
studied further. Parenthetically, we note that the locations of the nontrivial zeros of the 
Riemann zeta function appear to exhibit chaotic behaviour. If these positions are identified 


with energy level spacings, the corresponding potentials have many oscillations |14|. 


We would like to thank A. Gangopadhyaya and T. Imbo for valuable discussions and 
comments. U.S. acknowledges the hospitality of the Institute of Physics, Bhubaneswar 
where part of this work was done. This research was supported in part by the U.S. Depart¬ 
ment of Energy. 
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Figure Captions 


Figure 1: This diagram illustrates the energy levels of a cyclic shape invariant potential 
of period p. The lowest level is at zero energy, and the other excited energy levels are given 

by (B). 


/Q\ 

Figure 2: The superpotential Wq {x) and the corresponding symmetric potential 

/Q'\ 

V_ {x) for period p = 3 and a choice of energy level spacings loq = 1.0, ioi = 0.5, uj 2 = 2.3. 
The superpotential was obtained by a numerical solution of eqs. (O) by the Runge-Kutta 


method. The oscillation frequency given by eq. (^) is q^i = 1.097, and is in agreement with 
the figure. This is established by counting the number of maxima in any chosen interval of 
X values. The energy levels for this potential are given by expression (B) and are located at 
nujs, no's + 1.0, nojs + 1.5 with cus = 3.8 and n = 0,1, 2,... 


Figure 3: A plot of the superpotential for period p = 5 and energy level 

spacings loq = 0.8, cui = 1.5, U 2 = 1.0, cus = 3.0, UJ 4 = 0.2 . The asymptotic behavior 
has been subtracted in order to clearly reveal the two oscillation frequencies q^i = 0.4722 
and 952 = 2.001, predicted by eq. The maxima are spaced in agreement with the 

smaller frequency q^i whereas their “envelope” corresponds to the larger frequency 952 . 

Figure 4: This figure shows how the superpotential Wq^\x) — gets progressively 
more complicated and aquires more oscillation frequencies as the period p is successively 
doubled along the bifurcation route to chaos via the logistic map The parameter 

values chosen for the figures are: 

(a) r = 3.200 which gives period p = 2 and a stable 2-cycle with energy spacings 
ujo = 0.799, uji = 0.513 . 

(b) r = 3.500 which gives period p = A and a stable 4-cycle with ujq = 0.827, u)i = 

0.501, a ;2 = 0.875, ws = 0.383. The oscillation frequency from eq. (^) is qn = 0.323. 

(c) r = 3.561 which gives period p = 8 and a stable 8 -cycle with loq = 0.890, wi = 

0.348, 0 J 2 = 0.808, u!s = 0.552, UJ 4 = 0.881, UI 5 = 0.374, uq = 0.834, uj'i = 0.494. The three 

oscillation frequencies given by eq. are ggi = 0.1341, q ^2 = 0.3238, ggs = 0.7818. 
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